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Abstract 


The  compressible  hydrodynamic  approach  previously  developed  for  small  Mach 
number  non-cavitating  flows  has  been  extended  to  simulate  cavitating  flows  as  well  as 
non-cavitating  flows.  The  extension  is  made  possible  by  assuming  a  complex  equation  of 
state  relating  density  and  pressure  to  cover  the  liquid  phase  and  the  gas  phase.  Thus,  the 
cavitation  phenomenon  is  regarded  as  a  single-phase  flow  phenomenon  enabling  the 
elimination  of  the  cavity  closure  condition.  The  numerical  model  is  an  unsteady  3- 
dimensional  flow  model  based  on  a  large  eddy  simulation  approach.  It  is  applied  to  typical 
thin  hydrofoils  and  thick  hydrofoils  at  non-cavitating  conditions  and  various  cavitating 
flow  conditions,  including  moving  cavity,  stable  sheet  cavity  and  sheet  cavity/cloud  cavity 
cyclical  flow  conditions.  Computations  are  carried  out  primarily  for  2-dimensional  foils, 
but  3 -dimensional  flow  characteristics  are  also  examined. 

The  computational  results  are  compared  with  some  available  data;  good 
quantitative  and  qualitative  agreements  are  indicated.  It  is  considered  very  significant  that 
the  sheet  cavitation/cloud  cavitation  phenomenon  is  found  to  be  similar  to  the  viscous 
boundary  layer  flow  separation/vortex  shedding  and  washout  phenomenon  in  many 
respects.  Cavitation  is  found  to  trigger  boundary  layer  separation  in  otherwise  non- 
separated  flow. 
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I.  Introduction 


Traditionally,  cavitating  flows  are  regarded  as  two-phase  flows  and  fundamentally 
different  from  non-cavitating  flows.  Thus,  classical  works  of  Helmholtz  (1868)  and 
Kirchhoff  (1869)  neglect  the  density  of  the  gas  and  treat  cavitating  flows  as  free-surface 
potential  flows  of  the  liquid  phase.  This  tradition  has  been  extended  for  about  a  century 
by  numerous  investigators,  such  as  Riabouchinsky  (1919),  Wu  (1956),  and  Tulin  (1958), 
to  name  a  few.  The  most  serious  problem  limiting  the  progress  of  the  free-surface 
potential  flow  approach  is  the  need  for  a  cavity  closure  condition  for  finite  length  cavities 
to  approximate  the  real  fluid  flow  phenomenon.  Different  closure  models  have  been 
proposed  by  a  number  of  researchers,  such  as  the  image  model  (Riabouchinsky,  1919),  the 
re-entrant  jet  model,  the  transition  model  (Wu,  1956),  the  singularity  model  (Tulin,  1958), 
and  the  open  cavity  model  (Song,  1965).  The  free-surface  boundary  value  problem  can 
also  become  unmanageably  complex  when  there  are  many  cavitation  bubbles  or  when  a 
cavity  is  highly  unstable. 

To  overcome  the  difficulty  presented  by  the  unsteady  and  multiple  free  surfaces, 
Kubota  et  al.  (1992)  proposed  a  "bubble  two-phase  flow  model"  based  on  an  assumption 
that  the  fluid  consists  of  a  bubble-water  mixture;  and  the  bubble  volume  change  is 
governed  by  a  modified  Rayleigh's  equation.  In  this  way  the  existence  of  cavitation  is 
inferred  by  the  existence  of  a  large  void  fraction  region.  An  ingenious  idea— and  they  were 
able  to  qualitatively  mimic  the  appearance  of  unsteady  sheet  and  cloud  cavities.  However, 
the  paper  gives  no  quantitative  comparison  with  experimental  results.  Combining  the 
ideas  of  Kubota  et  al.  with  Song  and  Yuan's  (1988)  weakly  compressible  fluid,  a  simpler 
single-phase  cavitating  and  non-cavitating  flow  model  can  be  produced. 

The  weakly  compressible  flow  (also  known  as  compressible  hydrodynamic  flow) 
model  is  based  on  the  assumption  that  the  equation  for  state  of  water  can  be  approximated 
by  a  linear  equation  relating  pressure  and  density.  Recently  Song  (1996)  and  Song  and 
Chen  (1996)  showed  that  the  weakly  compressible  flow  equations  contain  the  inner 
solution  (compressibility  boundary  layer  flow)  and  the  outer  solution  (incompressible 
flow)  under  small  Mach  number  conditions.  Furthermore,  the  outer  solution  is 
independent  of  the  sound  speed  used  in  the  model;  therefore,  an  incompressible  flow 
solution  can  be  obtained  with  greatly  increased  speed.  This  approach  is  now  further 
extended  to  include  a  more  general  equation  of  state  that  approximates  the  pressure- 
density  relationship  extending  from  the  liquid  phase  to  the  gas  phase.  The  linear 
relationship  is  still  assumed  for  the  liquid  phase,  but  the  gas  phase  is  represented  by  a 
function  that  makes  the  density  decrease  very  rapidly  when  the  pressure  falls  below  a 
specified  critical  pressure.  The  fiilly  compressible  flow  equations  of  continuity  and 
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momentum  are  used  for  the  gas  phase,  but  all  quantities  are  continuous  at  the  gas-hquid 
interface. 

This  new  single  phase  unsteady  three-dimensional  flow  model  with  a  complex 
equation  of  state  is  applied  to:  1)  a  non-cavitating  flow  about  a  two-dimensional  foil,  2) 
an  E.N.  (Elliptic  Nose)  foil  with  sheet  and  cloud  cavitation,  and  3)  a  NACA0015  foil  with 
and  without  cavitation.  Quantitative  as  well  as  qualitative  comparisons  with  experimental 
data  are  made  and  good  agreements  are  obtained. 
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II.  Theoretical  Considerations 


The  importance  of  compressibility  has  long  been  overlooked  by  most 
hydrodynamicists.  But,  much  like  the  fact  that  viscosity  cannot  be  ignored  in  the  viscous 
(space)  boundary  layer.  Song  (1996)  has  shown  that  the  compressibility  cannot  be  ignored 
in  the  compressibility  (time)  boundary  layer.  The  fundamental  characteristics  of  the 
equations  of  motion  change  when  the  viscous  terms  are  dropped.  Similarly,  the 
characteristics  of  the  equation  of  continuity  changes  when  the  fluid  is  assumed  to  be 
incompressible  and  the  time  derivative  term  is  dropped.  Therein  lies  the  fundamental 
difficulty  in  representing  the  real  fluid  flows  with  the  classic  hydrodynamic  equations,  no 
matter  how  large  the  Reynolds  number  or  how  small  the  Mach  number.  The  remedy 
proposed  by  Song  and  Yuan  (1988)  is  to  assume  the  flow  to  be  "barotropic"  rather  than 
incompressible,  and  write  the  equation  of  state  as 


where  a  is  the  speed  of  sound. 


The  equation  of  continuity  can  now  be  written  as 


dt 


(1) 


(2) 


For  weakly  compressible  flow,  the  speed  of  sound  is  assumed  to  be  constant, 
hence,  the  density  is  a  linear  function  of  pressure.  In  this  case,  the  right-hand  side  of  Eq. 
(2)  is  zero  and  the  equation  of  continuity  simplifies  to 

^  +  (3) 

where  the  subscript  “o”  represents  a  reference  quantity.  Note  that  both  Eqs.  (2)  and  (3) 
are  written  in  conservative  form— one  with  and  the  other  without  a  source  term. 

For  a  cavitating  flow  problem,  it  is  necessary  to  simulate  at  least  approximately  the 
vaporization  phenomenon.  To  simulate  a  rapid  change  in  density  below  the  critical 
pressure,  the  following  equation  is  selected  for  the  equation  of  state. 
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for  Pe<P<Pc 


(4) 


p-  'Lap' 

i=0 

The  above  equation  is  joined  by  the  following  equation  at  the  upper  end, 

P-Po  =  i(p-Po),  for  Pc<P  (5) 

and  joined  by  the  following  equation  at  the  lower  end. 

P=Bp,  for  p<Pe  (6) 

The  coefficients  are  so  adjusted  that  the  resulting  pressure-density  curve  has  a  desirable 
shape. 


Since  the  Mach  number  of  the  flow  inside  a  cavity  may  not  be  small,  it  is  prudent 
to  use  the  equations  of  motion  for  fully  compressible  fluid  to  simulate  the  flow  inside  the 
cavity.  The  equations  of  continuity  and  motion  may  be  written  in  a  conservative  form  as 
follows. 


— +  V«F  =  5,  F  =  i£  +  JF  +  JtG  (7) 

dt 

where  S  is  source  which  may  or  may  not  exist  and  i,  j,  k  are  the  unit  vectors  in  x,  y,  and  z 
directions,  respectively.  The  functions  U,  E,  F,  G,  and  S  are  defined  in  the  following. 


u=\ 

T 

[p,pu,fyv,pwY 

(8) 

E  =  \ 

T 

\pa^u,pu^  +  p -  r ^,puv  -  T ^,puw  -  r 

(9) 

F  = 

rp 

[pa\puv-T^,pv^  +p-T^,pvw-T^^ 

(10) 

G  = 

T 

[yoa +  /? -  T„] 

(11) 

s-\ 

T 

yoV*  Va%0,0,o] 

(12) 
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III.  Numerical  Approach 


The  finite  volume  approach  with  MacCormack's  (1969)  predictor-corrector 
method  has  been  used  with  good  results.  Equation  (7)  is  first  integrated  over  a 
computational  finite  volume  and  divided  by  the  volume.  In  the  process,  the  volume 
integral  of  divergence  is  converted  to  the  surface  integral  by  applying  the  divergence 
theory.  The  resulting  equation  is 


dU  1  f  r  -JA  c 

—  + — \\AF»naA  =  S 
a  AV  ^ 


(13) 


In  the  above  equation  bar  means  volume  averaged  quantity,  A  is  surface  area,  and  n  is 
the  unit  normal  vector  on  the  surface  of  integration.  Since  the  quantities  to  be  calculated 
are  always  the  volume  averaged  values,  the  superscript  bar  will  be  dropped  hereafter.  The 
above  equation  is  a  simple  statement  of  the  law  of  conservation  where  the  rate  of  increase 

ofU  in  a  volume  is  equal  to  the  net  influx  of  F  plus  the  source  ^ .  Thus,  Eq.  (13)  can  be 
discretized  with  any  standard  finite  difference  scheme.  A  detailed  description  of  the 
process  used  herein  is  described  by  Yuan  (1988). 

Since  a  finite  volume  approach  can  resolve  only  the  quantity  of  scale  larger  than 
the  finite  volume  used,  it  is  necessary  to  use  a  closure  model  for  turbulent  flow  simulation. 
A  simple  Smagorinsky  (1963)  subgrid  scale  turbulence  model  is  used.  This  model 
assumes  that  the  shear  stress  tenser  %  can  be  calculated  by  adding  an  eddy  viscosity  Vt 
to  the  kinematic  viscosity  v.  Smagorinsky  proposed  the  following  equation  for  the  eddy 
viscosity. 


dXj  dXj 


(14) 


In  the  above  equation  C  is  a  coefficient  to  determined  and  A  is  the  size  of  the  finite 
volume.  Experience  indicates  that  a  C  distribution  which  is  zero  on  the  wall  and  increases 
to  0.12  outside  of  the  boundary  layer  generally  gives  satisfactory  results. 

The  size  of  the  finite  volume  used  is  often  too  large  to  adequately  resolve  the  sharp 
velocity  gradient  near  a  solid  wall.  For  this  reason  a  wall  function  or  a  partial  slip 
condition  (becomes  non-slip  when  boundary  layer  is  thick)  is  used;  see  He  and  Song 
(1991)  for  a  description  of  this  condition.  Because  the  cavity  is  assumed  to  be  a 
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continuous  extension  of  the  liquid,  no  free  surface  boundary  condition  is  necessary.  Only 
the  hydrodynamic  quantities  are  of  interest  (acoustics  are  ignored);  therefore,  an  arbitrarily 
large  value  of  M=0.4  is  chosen  for  the  liquid  part  of  the  flow  to  speed  up  the  computation. 
Recall  that  the  outer  solution  of  the  weakly  compressible  flow  equation  is  independent  of 
M  and  is  the  incompressible  flow. 
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IV.  Non-Cavitating  Flow  Over  A  Foil 


Flow  over  a  hydrofoil  even  at  non-cavitating  conditions  generates  a  very 
complicated  unsteady  flow  phenomena  similar  to  those  found  in  two-dimensional  diffuser 
flow  with  different  diffusion  angle,  especially  at  large  angle  of  attack.  At  a  small  angle  of 
attack,  the  flow  separation  may  occur  at  the  trailing  edge,  and  the  unsteady  characteristics 
of  the  flow  are  dominated  by  the  trailing  edge  vortex  shedding  much  like  the  unsteady 
shear-layer  flow.  As  the  attack  angle  increases,  the  unsteady  separation  point  on  the 
suction  side  moves  upstream.  At  a  large  enough  attack  angle,  the  separation  may  occur 
near  the  leading  edge,  and  a  periodical  vortex  shedding  like  the  well-known  Von  Karman's 
vortex  shedding  fi’om  a  bluff  body  may  take  place.  However,  in  some  cases,  a  number  of 
vortices  may  accumulate  near  the  foil  before  one  large  and  strong  negative  vortex  is 
washed  out.  Each  negative  vortex  being  washed  out  induces  a  positive  vortex  at  the 
trailing  edge  which  is  also  washed  out,  thus,  generating  a  low-ffequency  oscillation. 

The  vortex  shedding  and  vortex  washout  phenomena  described  above  are 
simulated  using  this  model,  and  the  results  are  presented  in  this  section.  The  main  effort 
of  the  non-cavitating  flow  study  is  focused  on  the  unsteadiness  of  the  flow  at  different 
attack  angles.  The  computed  mean  pressure  distribution  on  the  foil  agrees  very  well  with 
the  experimental  data.  Most  of  the  unsteady  characteristics  found  in  previous  experiments 
are  also  successfully  simulated. 


4-1  Previous  Studies 


Numerous  experiments  and  numerical  studies  of  the  flow  over  a  foil  at  non- 
cavitating  conditions  have  been  carried  out  for  several  decades,  and  a  good  understanding 
of  general  flow  patterns  has  been  achieved.  The  progress  and  the  state-of-the-art  on  the 
subject  can  be  learned  from  a  number  of  review  papers  (McCroskey,  1982;  Carr,  1988; 
and  Doligalski  et  al.,  1994).  However,  some  important  flow  characteristics,  such  as  the 
dynamics  of  stall  vortex  and  unsteady  flow  structures,  are  still  far  fi’om  being  completely 
understood.  Also  a  systematic  and  detailed  study  on  how  the  unsteady  stall  phenomena 
are  affected  by  the  attack  angle  change  is  lacking.  For  example,  Zaman  et  al.  (1989) 
observed  the  occurrence  of  low  fi-equency  oscillations  of  velocity  and  lift  when  the  flow 
approaches  the  stalling  condition.  But  some  questions  related  to  the  low  frequency 
oscillation  phenomenon  remain  unanswered  due  to  insufficient  experimental  data. 
Because  of  its  importance  to  rotor  craft,  low  frequency  stall  received  extensive 
experimental  studies  in  recent  years  by  Acharya  and  Metwall  (1992),  Currier  and  Fung 
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(1992),  and  Chandrasekhara  et  al.  (1993).  The  low  frequency  stall  frequently  that  occurs 
with  a  pitching  foil  has  also  been  found  to  occur  in  some  cases  around  a  stationary  foil. 

The  numerical  model  described  above  can  be  used  with  both  cavitating  and  non- 
cavitating  flow  cases.  In  fact,  when  the  cavitation  number  is  set  to  a  very  large  number, 
the  flow  becomes  non-cavitating  flow.  In  this  section,  the  numerical  model  is  used  to 
simulate  the  flow  over  a  stationary  foil  at  different  attack  angles  under  non-cavitating  flow 
conditions.  The  objectives  of  the  non-cavitating  flow  study  are:  1)  demonstrating  how  the 
unsteadiness  of  the  vortex  generation  process  is  affected  by  the  attack  angle;  2)  revealing 
more  details  of  the  dynamic  stall  process  by  means  of  numerical  visualization;  3)  providing 
the  bases  for  comparing  and  facilitating  a  fiirther  understanding  of  cavitating  flow  around 
hydrofoils. 

4-2  Model  Validation 

Very  small  Mach  number  flow  or  a  weakly  compressible  flow  is  considered,  and  p 
and  a  in  Eqs.  (8-12)  are  set  equal  to  constants.  A  body-fitted  grid  system  is  used  to  ensure 
accuracy.  To  assure  that  the  computed  results  will  be  grid  independent,  three  mesh 
systems  (191x32x3,  381x64x3,  761x128x3)  were  tested.  This  model  contains  only  three 
grid  points  in  the  spanwise  direction  so  that  basically  two-dimensional  flow  is  represented. 
Part  of  the  medium  mesh  system  (381x64x3)  around  the  NACA0015  foil  is  shown  in  Fig. 
4-1.  It  was  found  that  the  medium  mesh  system  is  fine  enough  to  resolve  the  large-scale 
unsteady  structures.  Therefore,  the  results  discussed  in  this  report  are  based  on  medium 
mesh  for  comparison  purposes  unless  otherwise  indicated.  The  Reynolds  number  used  in 

the  compution  is  10^- 

The  numerical  model  was  first  tested  on  NACA0015  standard  foil,  which  has  been 
intensively  studied  by  Izumida  et  al.  (1980)  and  Kubota  et  al.  (1992)  by  experiment  and 
computation  under  cavitating  and  non-cavitating  conditions.  Fig.  4-2  shows  the 
comparison  of  the  computed  time-averaged  pressure  distribution  on  the  foil  surface  with 
the  experimental  data  of  Izumida  et  al.  for  attack  angle  (a)  of  0,  4,  and  8  degrees.  The 
experimental  data  were  taken  only  on  the  suction  side.  As  the  figure  indicates,  the 
computed  results  agree  very  well  with  the  experiments.  The  computed  lift  coefficients  at 
three  attack  angles  are  0.0,  0.474,  and  0.905,  respectively. 

To  test  the  effect  of  the  grid  size,  three  mesh  systems  (191x32x3,  381x64x3,  and 
761x128x3)  were  used  for  the  case  of  a  =8°.  The  time-averaged  pressure  distributions 
based  on  the  three  meshes  are  shown  in  Fig.  4-3.  It  appears  that  the  results  using  the 
medium  size  mesh  are  very  close  to  those  using  the  fine  mesh.  Since  the  fine  mesh  takes 
more  than  eight  times  the  work  to  generate  the  same  results  than  the  medium  mesh 
system,  the  computed  results  discussed  below  are  all  based  on  the  medium  mesh  system. 
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4-3  Unsteady  Structures 

A  non-symmetrical  foil  NACA662-415  was  chosen  for  the  unsteady  flow  study. 
Four  different  attack  angles  6°,  10°  ,  15°  ,  and  20°  were  selected  for  the  study.  Results 
showed  that  the  flows  at  these  attack  angles  represent  different  unsteady  flow  regimes. 
These  flow  regimes  are  very  similar  to  the  flow  regimes  in  two-dimensional  S5mimetrical 
diflusers  of  different  diffusion  angles,  as  found  experimentally  by  Fox  and  Kline  (1962) 
and  computationally  by  He  and  Song  (1991). 

4.3.1  Trailing  edge  unsteadiness  ( a=6 °) 

At  this  small  attack  angle  the  flow  separates  near  the  trailing  edge  on  the  suction 
side,  which  results  in  unsteady  vortex  shedding  from  the  trailing  edge.  The  separation 
point  moves  upstream  as  the  attack  angle  increases.  Fig.  4-4  shows  instantaneous 
vorticity  contours  during  one  period  of  the  trailing-edge  vortex  shedding.  In  the  other 
region  away  from  the  trailing  edge,  the  flow  is  steady,  and  there  is  no  leading  edge 
separation. 

Fig.  4-5  shows  the  power  spectra  of  lift  and  velocity  component  (u)  recorded  at 
five  positions,  as  sketched  at  the  top  of  the  figure.  As  the  figure  indicates,  only  the  lift  and 
the  velocity  in  the  wake  near  the  trailing  edge  (position  4)  show  the  unsteadiness  due  to 
the  trailing-edge  vortex  shedding.  The  Strouhal  number  (fc-sin  a/U)  of  the  wake 
oscillation  is  0.60.  This  flow  regime  corresponds  to  the  steady  flow  regime  of  the  two- 
dimensional  diffuser  at  a  6°  diffusion  angle  (see  He  and  Song,  1991). 

4.3.2  Leading  edge  separation  (a=10°) 

When  the  attack  angle  is  about  10°,  the  leading  edge  separation  occurs,  which 
generates  a  leading  edge  separation  bubble  (clockwise  rotating  eddy).  Vorticity  generated 
at  the  separation  point  continuously  feeds  the  eddy,  causing  the  eddy  to  grow  in  size  and 
intensity.  When  the  eddy  grows  to  a  certain  size,  part  detaches  and  forms  a  primary 
dynamic  stall  vortex.  Such  a  process  of  vortex  generation  and  detachment  is  similar  to 
bluff-body  vortex  shedding.  The  detached  vortex  moves  downstream  with  the  main  flow 
(U)  at  an  approximate  speed  (Doligalski  et  al.,  1994), 

V=U-—  (15) 

2d 

where  k  is  the  vortex  strength,  and  d  is  the  distance  from  the  center  of  the  vortex  to  the 
foil  surface.  Due  to  the  relatively  small  vortex  strength,  the  speed  of  vortex  convection  is 
not  much  different  from  that  of  the  main  flow. 

Fig.  4-6  shows  the  instantaneous  vorticity  contours  during  one  period  of  the 
vortex  shedding  at  a=10°.  The  detached  vortex  remains  very  close  to  the  foil  surface 
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while  being  transported  from  the  leading  edge  to  the  trailing  edge.  When  this  clockwise 
rotating  vortex  reaches  the  trailing  edge,  it  induces  a  counterclockwise  rotating  vortex  at 
the  trailing  edge.  Four  instantaneous  close-up  views  of  the  leading  edge  vortex  within  one 
cycle  of  shedding  is  shown  in  Fig.  4-7. 

Fig.  4-8  shows  the  power  spectra  of  lift  and  the  velocity  component  (u)  recorded 
at  the  same  five  positions,  as  stated  above.  The  most  significant  spectral  peak  in  the  figure 
is  at  the  leading  edge,  where  the  unsteadiness  due  to  the  vortex  shedding  is  very  strong. 
The  Strouhal  number  (fc-sin  oAJ)  of  the  leading  edge  vortex  shedding  is  0.2.  This  flow 
regime  corresponds  to  the  small  transitory  stall  regime  of  two-dimensional  diffuser  at  a 
13°  diffusion  angle. 

4.3.3  Large  vortex  washout  (a= 15°) 

As  the  attack  angle  increases,  the  shed  vortex  becomes  stronger,  and  its  moving 
speed  decreases  according  to  Eq.  (15).  At  this  large  attack  angle,  a  vortex  merging  and 
washout  phenomenon  take  place.  A  vortex  that  was  shed  earlier  moves  slowly  and 
appears  to  wait  for  the  trailing  vortices  to  catch-up  and  merge,  resulting  in  a  large  vortex 
in  the  vicinity  of  the  foil.  This  large  combined  vortex  is  eventually  washed  out  of  the  foil 
region  into  the  wake.  The  large  stall  washout  generates  strong,  low  frequency  oscillatory 
lift  and  drag. 

Fig.  4-9  shows  one  period  of  combined  vortex  washout  featured  by  the 
instantaneous  vorticity  contours  at  a=15°.  In  Fig.  4-9(a),  one  large  stall  vortex  is  washed 
out  from  the  foil.  The  vortex  at  the  trailing  edge  in  this  figure  is  induced  by  the  large 
clockwise  stall  vortex.  Hence,  the  induced  small  but  high  strength  vortex  at  the  trailing 
edge  at  this  stage  is  a  counterclockwise  vortex.  In  the  leading  edge  region,  as  more  shed 
vortices  shed  and  merge,  the  combined  vortex  becomes  larger  and  larger.  Since  the  large 
vortex  moves  downstream  slowly,  the  new  shed  vortices  from  the  leading  edge  can  still 
catch  up  before  it  reaches  the  trailing  edge.  This  flow  regime  corresponds  to  the  large 
transitory  stall  of  the  two-dimensional  diffuser  of  18°  diffusion  angle  (He  and  Song, 
1991).  Such  a  vortex  structure  and  its  unsteady  behavior  are  also  very  similar  to  the 
coalescence  of  vortex  cavitation  observed  in  a  water  tunnel  experiment  using  high-speed 
photograph  (Kubota,  et  al.,  1992).  In  fact,  a  recent  study  (Song,  et  al.,  1997)  to  be 
described  later  shows  that  the  cavitation-induced  unsteady  flow  structure  over  a  foil  at  a 
small  attack  angle  is  similar  to  that  of  non-cavity  flow  at  a  larger  angle. 

At  the  same  five  positions  as  before,  the  instantaneous  velocity  component  (u)  and 
lift  were  recorded.  Their  power  spectra  are  shown  in  Fig.  4-10.  Due  to  the  larger  leading 
edge  separation  eddy.  Position  1  falls  inside  of  the  eddy  and  has  very  little  velocity 
fluctuation.  In  contrast.  Position  2  is  close  to  the  end  of  the  separation  eddy  and  has 
strong  oscillation  of  the  leading  eddy  shedding  frequency  (fi*c-sin  aAJ=0. 17).  At  Position 
5  (near  wake  region)  the  low  frequency  large  eddy  washout  (f2’C*sin  a/U=0.05)  is  more 
significant.  An  interesting  feature  is  that  the  unsteady  lift  contains  two  dominant 
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frequencies:  the  low  washout  frequency  and  a  higher  frequency  (fs  -c-sin  (x/U-0.3).  The 
higher  frequency  unsteadiness  is  found  to  be  related  to  the  trailing  edge  secondary  vortex 
shedding  similar  to  that  at  low  angle  of  attack.  The  unsteadiness  of  the  leading  edge 
vortex  shedding  does  not  contribute  much  to  lift  oscillation  because  there  is  more  than 
one  such  vortex  on  the  foil,  and  their  effect  on  lift  cancel  out  each  other. 

Fig.  4-11  shows  the  details  of  the  leading  edge  vortex  shedding  featured  by  the 
velocity  vectors.  The  unstable  large  recirculation  at  the  leading  edge  is  clearly  seen.  The 
complex  flow  structure  contains  the  spike-like  eruption,  shear  layer  roll-up,  and  viscous- 
inviscid  interaction,  as  described  by  Dohgalski  et  al.  (1994),  and  visualized  in  experiments 
by  Acharya  and  Metwally  (1992).  A  similar  low  frequency,  combined  eddy  washout 
phenomenon  was  also  found  in  an  experiment  by  Zaman  et  al.  (1989).  The  shedding 
vortex  accumulation  on  the  foil  surface  is  also  very  similar  to  the  cavitation  coalescence 
found  in  experiments  using  high-speed  photographs  (Kubota  et  al.,  1992). 

4.3.4  High  attack  angle  stall  (a=20°) 

As  the  attack  angle  is  increased  further,  the  shed  vortex  from  the  leading  edge  is 
more  detached  from  the  foil  surface.  As  a  result,  the  vortex  merging  becomes  less 
significant.  The  extreme  case  is  a=90°,  in  which  case  the  only  large-scale  unsteadiness  is 
the  bluff-body  vortex  shedding. 

Fig.  4-12  shows  one  period  of  the  leading-edge  vortex  shedding  featured  by  the 
instantaneous  vorticity  contours  at  a=20°.  The  power  spectra  of  u  at  the  five  positions 
and  the  lift  are  displayed  in  Fig.  4-13.  It  is  seen  that  the  unsteady  flow  velocity  at  all  five 
positions  are  dominated  by  the  leading-edge  vortex  shedding,  which  has  a  Strouhal 
Number  (fc-sin  oAJ)  of  0.17,  although  large  eddy  washout  can  still  occur  from  time  to 
time.  In  this  case,  the  lift  power  spectrum  shows  a  small  peak  representing  the  leading- 
edge  vortex  shedding. 

In  summary,  the  above  discussed  four  cases  represent  four  flow  regimes  of  flow 
over  a  hydrofoil  at  different  attack  angles.  Each  regime  features  different  large-scale 
unsteady  structures.  The  flow  regimes  are  very  similar  to  those  found  in  a  two- 
dimensional  diffuser  flow  with  different  opening  angle  (Fox  and  Kine,  1962;  He  and  Song, 
1991). 
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V.  Cavitating  Flow  Over  A  Thin  Foil 


As  mentioned  before,  flow  at  a  cavitating  condition  is  much  more  complicated  than 
the  corresponding  non-cavitating  flow.  Not  only  do  cavity  generation  and  collapse  lend 
difficulties  for  numerical  simulation,  but  there  are  more  factors  to  be  considered.  In  the 
current  cavitating  flow  study,  two  types  of  foil  shapes  were  examined:  the  NACA0015 
which  represents  a  t5T)ical  thick  foil  (thickness  ratio  of  12%);  and  the  E.N.  (Elliptic  Nose) 
foil  which  is  a  relatively  thin  foil.  The  cavitating  flow  results  over  the  NACAOOl  5  foil  will 
be  discussed  separately  in  the  next  section  since  different  phenomena  were  observed.  In 
this  section,  the  cavitating  flow  results  over  an  E.N.  foil  is  analyzed. 

It  is  known  that  during  non-cavitating  conditions,  the  separation  point  at  the 
suction  side  of  a  thick  foil  moves  upstream  to  the  half  foil  chord  as  the  attack  angle 
increases  before  the  leading  edge  vortex  shedding  occurs.  However,  for  a  thin  foil,  the 
leading  edge  vortex  shedding  occurs  at  a  relatively  small  attack  angle  due  to  the  high  peak 
negative  pressure  at  the  leading  edge.  The  so-called  E.N.  foil  is  a  symmetrical  hydrofoil 
with  a  thickness  ratio  of  8%.  The  foil  geometry  can  be  expressed  as  (Kawanami  et  al., 
1996), 


= +(0.05746  x- 0.06079) 

(0.4187  <x  <  l.O) 

(16) 

.y  =  ±:^A/x(0.6-x) 

VI 

X! 

VI 

(17) 

The  other  reason  to  study  this  foil  is  that  an  extensive  experimental  study  on  cloud 
cavitation  around  this  foil  and  its  unsteady  structures  were  reported  by  Kubota  et  al. 
(1989).  Recently,  Kawanami  et  al.  (1996)  reported  a  follow-up  study  of  this  foil  about  the 
mechanism  and  control  of  cloud  cavitation.  In  order  to  compare  the  modeling  results, 
conditions  were  set  to  be  the  same  as  those  used  in  their  experiments. 

The  same  mesh  system  used  for  the  non-cavitating  flow  computation  was  used. 
To  compare  the  experiments,  an  attack  angle  of  6.2  degrees  was  used  in  the  computation, 
but  different  cavitation  numbers  were  tested. 


5.1  Non-Cavitating  Flow  Results 

As  the  basis  for  comparison  with  cavitating  flow  characteristics,  non-cavitating 
flow  at  the  same  modeling  configuration  was  first  simulated.  It  was  found  that  the  flow 
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field  over  the  E.N.  foil  at  the  attack  angle  of  6.2  degrees  at  non-cavitating  conditions  is 
mostly  steady  except  for  the  unsteady  vortex  shedding  at  the  trailing  edge. 

Fig.  5-1  shows  one  period  of  the  vortex  shedding  at  the  trailing  edge  featured  by 
the  vorticity  contours.  As  indicated  in  the  figure,  the  shed  vortices  from  the  trailing  edge 
gradually  defuse  in  the  wake  region.  There  is  no  separation  at  the  leading  edge.  Fig.  5-2 
shows  the  time-averaged  velocity  vectors.  No  boundary  layer  separation  except  at  the 
trailing  edge  is  shown.  The  close-up  of  the  velocity  vectors  at  the  leading  edge  is  shown 
in  Fig.  5-3.  The  time-averaged  pressure  coefficient  distribution  around  the  foil  is  shown  in 
Fig.  5-4.  Fig.  5-5  displays  the  pressure  coefficient  distribution  on  the  foil  surface.  At  the 
leading  edge,  there  is  a  sharp  negative  pressure  peak  higher  than  4.  The  sharp  negative 
peak  is  typical  in  a  thin  foil  like  this.  The  computed  time-averaged  lift  coefficient  at  this 
attack  angle  (6.2°)  is  0.754. 


5.2  Cavitating  Flow  Results 

At  cavitating  conditions  under  the  same  attack  angle  (a=6.2°)  the  flow  patterns  are 
totally  different  fi'om  those  of  the  non-cavitating  flow.  When  the  ambient  pressure  is 
lowered  below  the  critical  value,  cavitation  starts  in  the  leading  edge  region  due  to  the 
peak  negative  pressure  at  the  leading  edge.  The  cavitation  bubble  exhibits  a  distinctively 
periodical  behavior,  alternating  between  the  growth  stage  and  the  break-up  stage.  During 
the  growth  stage,  vorticity  is  continuously  produced  near  the  leading  edge  and  fed  into  the 
bubble  making  the  bubble  grow  and  form  a  large  sheet  cavity.  The  sheet  cavity  is 
relatively  stable  during  the  growth  stage.  But  it  becomes  unstable  and  breaks  away  fi’om 
the  upper  section  and  forms  a  large  clockwise  rotating  free  eddy.  Actually,  this  fi’ee  eddy 
is  highly  turbulent  and  contains  numerous  small  cavitating  eddies  called  “cloud  cavitation.” 

Fig.  5-6  shows  one  cycle  of  the  cloud  cavitation  featured  by  instantaneous  vorticity 
contours  when  the  cavitation  number  is  about  1.2.  In  Fig.  5-6(a),  a  large  sheet  cavity  is 
still  growing  as  one  cloud  cavitation  has  just  been  washed  out  away  fi’om  the  foil’s  trailing 
edge.  In  Fig.  5-6(b),  a  large  portion  of  the  sheet  cavity  is  beginning  to  break  away  at 
about  the  1/4  chord  point.  In  Fig.  5-6  (c)  the  sheet  cavity  that  has  broken  away  turned 
into  a  full  fledged  cloud  cavity  and  moved  to  about  the  3/4  chord  point  where  relatively 
high  ambient  pressure  is  causing  the  cloud  cavity  to  start  collapsing.  From  Fig.  5-6(d)  to 
Fig.  5-6(e),  the  collapsed  cloud  cavitation  is  gradually  dying  out  as  it  is  washed  away  from 
the  foil  section,  while  another  sheet  ca'vity  is  forming.  This  completes  one  cycle  of  sheet 
cavity  growth-break-up  into  cloud  cavity— collapse  and  washout  of  cloud  cavity.  This 
same  periodical  characteristics  of  sheet  and  cloud  cavitation  was  observed  in  experiments 
by  Kubota  et  al.  (1989).  Fig.  5-7  shows  the  photographs  of  one  period  of  the  cavitation 
cloud  fi’om  generation  to  washout,  taken  from  the  experiment  by  Kubota  et  al.  (1989). 
Although  Fig.  5-6  and  Fig.  5-7  show  flow  visualization  using  different  physical  quantities, 
the  flow  patterns  shown  are  strikingly  similar.  The  computed  velocity  field  also  illustrates 
the  same  unsteady  characteristics  as  shown  in  Fig.  5-8(a).  Similar  velocity  measurements 
during  one  cycle  of  cloud  cavitation  based  on  LDV  measurements  of  Kubota  et  al.  are 
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shown  in  Fig.  5-8(b).  It  appears  that  the  large-scale  unsteady  structure  of  the  cloud 
cavitation  obtained  by  the  computational  model  agrees  very  well  with  the  experimental 
results.  Compared  with  the  non-cavitating  flow  shown  in  Fig.  5-1,  cavitating  flow  is  much 
more  unstable.  But  there  are  many  similarities  when  compared  with  the  non-cavitating 
flow  at  an  attack  angle  of  15°  where  vortex  shedding  and  washout  occur  as  shown  in  Fig. 
4-9. 


The  velocity  vectors  around  the  foil  at  the  moment  when  the  sheet  cavity  begins  to 
break  away  from  the  foil’s  upper  section,  corresponding  to  Figs.  5-6(c),  is  shown  in  Figs. 
5-9.  This  figure  indicates  the  three  different  vortex  structures  as  mentioned  above.  They 
are  the  leading  edge  vortex  production  which  takes  place  in  the  first  10%  of  the  foil 
section  (between  points  1  and  2),  the  large  sheet  cavity  which  extends  to  25%  of  the  foil 
(between  points  2  and  3),  and  the  cavitation  cloud  which  moves  from  the  foil’s  middle 
section  to  the  wake  (downstream  of  point  3).  The  highly  unsteady  leading  edge  vortex 
production  zone  should  probably  be  regarded  as  a  part  of  the  sheet  cavity  structure.  Fig. 
5-10  shows  the  close-up  of  the  velocity  field  in  the  vorticity  production  zone  (between 
points  1  and  2)  during  one  cycle  of  oscillation. 

Time-space  (t,y),  distributions  of  vorticity  at  three  sections  (x/c=0.25,  0.65,  and 
1.05)  are  shown  in  Fig.  5-11.  At  x/c=0.25.  Fig.  5-1 1(a),  the  vorticity  contour  lines  are 
nearly  parallel  to  the  t-axis,  indicating  that  the  sheet  cavity  is  relatively  stable.  Periodically 
near  the  time  of  break  up,  the  contour  lines  become  normal  to  the  time  axis  indicating 
rapid  change.  At  x/c=0.65  as  shown  in  Fig.  5-1 1(b),  the  large  vortex  structure  of  cloud 
Cavitation  appears  periodically.  After  each  cloud  cavitation  vwth  negative  vorticity  passing 
this  foil  section,  a  small  positive  vortex  structure  is  induced.  In  Fig.  5-1 1(c),  located  0.05c 
downstream  of  the  trailing  edge,  the  vortex  of  cloud  cavitation  becomes  weaker.  Positive 
vorticity  induced  by  the  flow  on  the  pressure  side  is  present  near  y  =  0  most  of  the  time. 
Therefore,  the  vorticity  distribution  at  this  wake  section  changes  from  negative  in  the 
upper  level  to  positive  in  the  lower  level.  The  computed  vortex  structure  at  the  wake 
section  was  also  recorded  in  Kubota  et  al.'s  experiment  (1989)  using  LDV.  Fig.  5-12 
shows  the  experimental  results. 

Fig.  5-13  (a)  shows  the  computed  pressure  fluctuations  recorded  at  three  locations 
(x/c=0.45,  0.65,  and  0.85),  respectively.  At  x/c=0.45  as  shown  in  Fig.  5-7(a),  the  high 
frequency  pressure  oscillations  with  sharp  cut-off  at  p  =  0  represents  the  sheet  cavity 
pressure  fluctuations.  At  x/c=0.85  as  shown  at  the  top  of  Fig.  5- 13  (a),  the  pressure 
fluctuation  reflects  mostly  the  change  due  to  the  cloud  cavitation.  The  computed 
frequency  is  about  23  Hz  and  the  maximum  pressure  coefficient  is  about  0.58xl0E5  Pa. 
The  measured  pressure  data  at  the  same  locations  (Kubota  et  al.,  1989)  are  shown  in  Fig. 
5- 13(b).  The  measured  frequency  is  25.4  Hz  and  the  maximum  pressure  coefficient  is 
about  0.52xl0E5  Pa.  The  agreement  between  the  computation  and  the  experiment  is 
surprisingly  good.  For  a  longer  time  period,  the  pressure  change  with  time  at  five 
locations  along  the  foil  surface  is  shown  in  Fig.  5-14.  At  x/c=0.25,  high  frequency 
oscillation  of  the  sheet  cavity  is  predominant,  but  the  dominant  frequency  gradually 
decreases  as  x/c  increases.  At  x/c=1.05,  the  pressure  fluctuation  is  fully  dominated  by  the 
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large-scale  cavitation  cloud.  The  computed  normal  velocity  at  x/c=0.25  near  the  foil 
surface  (y/c=0. 133)  is  shown  in  Fig.  5-15(a).  The  corresponding  data  of  Kubota  et  al.  are 
shown  in  Fig.  5- 1 5(b). 

Fig.  5-16  shows  the  computed  time-averaged  pressure  distribution  on  the  foil 
surface  at  different  cavitation  numbers  compared  with  that  of  non-cavity  flow.  As  the 
cavitation  number  decreases,  the  flat  pressure  zone  on  the  suction  side,  which  represents 
the  vapor  pressure,  becomes  longer,  indicating  that  the  cavity  is  getting  longer. 
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VI.  Cavitating  Flow  Over  A  Thick  Foil 


As  noted  before,  the  cavitation  characteristics  of  a  thick  foil  are  somewhat 
different  from  those  of  a  thin  foil.  The  location  of  cavitation  occurrence  varies  more  with 
the  attack  angle  and  the  phenomenon  of  sheet  cavity  breakdown,  leading  to  less 
predictability  of  cloud  cavitation.  In  this  section,  cavitating  flow  results  over  a  typical 
thick  foil  NACA0015  are  presented.  The  same  mesh  system  described  previously  was  also 
used  in  this  case  to  study  the  basic  characteristics  of  the  cavitating  flow.  In  addition,  the 
number  of  grids  in  the  spanwise  direction  has  been  increased  from  3  to  15  and  then  to  30 
to  study  the  three-dimensionality  of  the  flow. 

As  shown  before,  non-cavitating  flows  over  a  NACA0015  foil  with  various  attack 
angle  were  simulated  first  to  check  for  the  accuracy  of  the  model  as  well  as  to  serve  as  a 
reference.  The  comparison  of  computed  pressure  coefficient  distribution  on  the  foil 
surface  and  the  experimental  data  was  shown  in  Fig.  4-2.  Good  agreements  have  been 
noted.  Fig.  6-1  shows  the  pressure-density  relationship  in  the  gas  phase.  Note  that  the 
critical  pressure  coefficient  of  0.10  can  be  chosen  arbitrarily.  This  figure  shows  a  rapid 
drop  of  density  to  10  %  of  water  within  the  pressure  coefficient  drop  of  0.01,  which  is 
substantially  milder  than  the  real  vaporization  process.  But  the  numerical  results  of 
cavitation  phenomenon  appear  insensitive  to  this  curve.  The  equation  of  state  for  Cp  > 
0. 10  is  assumed  to  be  linear  and  matches  this  curve  at  the  critical  point. 


6.1  Comparison  with  Experiments 

Only  a  high  speed  photograph  of  cavitating  flow  at  a  20  degree  attack  angle  and 
the  cavitation  number  of  1.7  published  by  Kubota  et  al.  (1992)  is  available  for 
comparison  for  this  foil.  Fig.  6-2  (a)  is  the  photograph  of  the  cavitating  flow  in  a  water 
tunnel,  and  Fig.  6-2  (b)  is  the  corresponding  computed  instantaneous  density  distribution. 
Here  the  cavitation  zones  indicated  by  the  relative  density  less  than  one  are  painted  white 
and  the  non-cavitating  zones  are  shown  in  black.  Small  amounts  of  sheet  cavitation  near 
the  nose  and  a  larger  cloud  cavitation  away  from  the  trailing  edge  are  captured  both 
experimentally  and  computationally;  they  show  striking  similarities. 
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6.2  Qualitative  Description  of  Two-Dimensional  Cavity  Flows 

6. 2. 1  Moving  and  fixed  (sheet)  cavitation  at  small  attack  angle 

As  discussed  in  Sections  V  and  VI,  two-dimensional  flows  are  actually  computed 
with  the  three-dimensional  model  but  using  only  three  grids  in  the  spanwise  direction. 
Conditions  at  various  attack  angles  and  cavitation  numbers  were  simulated  and  the 
pressure,  density,  velocity,  and  vorticity  output  were  studied.  Computed  instantaneous 
pressure  fields  for  the  case  of  a=4°  and  a  =1.5  are  depicted  in  Fig.  6-3.  Corresponding 
density  and  vorticity  distributions  are  shown  in  Fig.  6-4  and  Fig.  6-5,  respectively.  The 
first  two  figures  indicate  the  existence  of  a  moving  cavity  that  first  appears  near  the 
midchord  that  expands  and  contracts  as  it  moves  downstream  along  the  foil  surface. 
When  all  three  figures  are  compared,  it  is  clear  that  the  cavitation  bubble  originates  where 
the  viscous  boundary  layer  separation  first  takes  place  and  then  coincides  with  the 
trajectory  of  the  center  of  the  breakaway  eddy.  Since  there  would  be  no  boundary  layer 
separation  under  non-cavitating  conditions  at  the  same  attack  angle,  it  is  also  clear  that,  in 
this  case,  cavitation  induces  boundary  layer  separation  rather  than  the  other  way  around. 

The  simulated  pressure,  density,  and  vorticity  distributions  at  the  same  attack  angle 
but  at  smaller  cavitation  number  (o  =1.0)  are  shown  in  Figs.  6-6,  6-7,  and  6-8, 
respectively.  A  major  change  in  flow  pattern  has  taken  place  from  the  previous  case.  A 
fixed  cavity  starting  near  the  nose  and  extending  to  about  1/3  to  1/2  chord  position  is 
apparently  taking  place.  One  vortex  sheet  starting  from  the  nose  and  another  sheet 
starting  from  the  tail  of  the  cavity  can  be  observed  from  Fig.  6-8.  The  last  picture  in  (a) 
and  (b)  appears  to  show  that  the  negative  pressure  wave  may  propagate  around  the 
trailing  edge  and  cause  temporary  cavitation  on  the  pressure  side  of  the  foil.  This 
phenomenon  may  not  be  observable  in  a  physical  experiment  because  the  actual  sound 
speed  is  about  100  times  as  fast  as  the  value  used  for  simulation. 

6. 2. 2  Intermittent  occurrence  of  cloud  cavitation  at  medium  attack  angle 

At  attack  angle  of  8  degrees,  the  sheet  cavity  produced  at  small  cavitation  number 
is  quite  unstable  and  may  break  down  to  produce  cloud  cavitation.  But  unlike  the  case  of 
EN  foil  described  before,  the  occurrence  of  cloud  cavitation  is  not  regular  for  NACA0015 
foil.  Figs.  6-9,  6-10,  and  6-11,  respectively,  show  the  pressure,  density,  and  vorticity 
fields  at  o  =1.0  in  a  time  segment  when  sheet  cavity  is  not  producing  cloud  cavity.  The 
flow  patterns  shown  in  these  figures  are  very  similar  to  the  flow  patterns  shown  in  Figs. 
6-6,  6-7,  and  6-8  which  are  for  the  case  of  4  degree  attack  angle  but  at  the  same  cavitation 
number.  For  exactly  the  same  attack  angle  and  cavitation  number  but  at  different  time 
segment,  the  sheet  cavity  can  become  unstable  and  produces  cloud  cavity  as  indicated  by 
Figs.  6-12,  6-13,  and  6-14.  The  occurrence  of  cloud  cavity  appears  to  be  rather  random 
and  not  to  follow  any  deterministic  rule. 

When  we  observe  the  change  in  pressure  field  with  time  as  shown  in  Fig.  6-12 
from  (a)  to  (f),  it  is  possible  to  see  the  collapsing  cloud  cavity  generating  high  pressure 
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concentration  near  the  1/3  chord  position  as  indicated  by  a  dark  area.  This  pressure  is  far 
above  the  stagnation  pressure,  indicating  the  dynamic  and  compressibility  effect  of  a 
collapsing  cavity.  From  (d)  to  (f),  it  is  possible  to  observe  the  propagation  of  the 
compression  wave  away  from  the  center  of  the  collapse.  This  pressure  wave  is  not 
quantitatively  accurate  because  the  sound  speed  used  in  the  computation  is  only  about 
100th  of  the  real  sound  speed. 

6. 2. 3  Continuous  occurrence  of  cloud  cavitation  at  large  attack  angle 

At  large  attack  angle  and  low  cavitation  number,  even  a  thick  foil  as  the  NACA 
0015  will  continuously  produce  cloud  cavitation  or  supercavitation.  The  instantaneous 
patterns  of  pressure,  density  and  vorticity  for  the  case  of  p=20°  and  a  ==1.0  are  shown  in 
Figs.  6-15,  6-16,  and  6-17,  respectively. 

6.3  Three  Dimensionality  of  Cavitating  Flow 

To  study  the  three  dimensionality  of  the  flow  the  finite  span  equal  to  a  0.7  chord 
length  was  selected  and  the  span  was  divided  into  30  equal  grid  intervals.  A  nonslip 
boundary  condition  was  used  to  simulate  a  typical  water  tunnel  testing  condition.  Only 
the  results  for  the  case  of  an  8°  attack  angle  and  a  cavitation  number  of  1.0  at  an  instant 
when  cloud  cavity  exists  are  presented  here. 

Fig.  6-18  shows  the  pressure  distribution  over  three  different  cross-sections:  one 
near  the  left  wall,  one  at  the  center  of  the  span,  and  the  other  near  the  right  wall.  The 
pressure  patterns  are  very  different,  and  the  collapsing  cloud  cavity  condition  can  be  seen 
only  near  the  right  wall  at  this  instant.  The  corresponding  density  distribution  is  shown  in 
Fig.  6-19;  the  velocity  distribution  is  shown  in  Fig.  6-20. 

Fig.  6-21  shows  the  pressure  distributions  over  three  surfaces  roughly  parallel  with 
the  foil  suction  surface  but  at  different  distances  away  from  the  foil.  This  figure  again 
clearly  shows  that  the  cavity  collapse  occurs  near  the  right  wall  and  near  the  foil  surface. 
The  corresponding  density  distribution  and  velocity  distribution  are  shown  in  Figs.  6-22 
and  6-23 .  Note  that  the  upstream  half  of  Section  1  in  Fig.  6-23  is  the  foil  surface  and  the 
velocity  is  zero  according  to  the  boundary  condition.  A  wavy  vortex  pattern  in  the  span 
direction  is  a  typical  flow  pattern  over  a  finite  span  foil  with  or  without  cavitation. 
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VII.  Conclusions 


The  Compressible  Hydrodynamic  Approach  previously  developed  for  large 
Reynolds  number  -  small  Mach  number  flows  has  been  extended  and  applied  to  cavitating 
and  non-cavitating  flows  about  hydrofoils  of  different  shapes  at  different  attack  angles  and 
for  different  cavitation  numbers.  The  unifying  factor  between  the  previous  approach  and 
the  extended  approach  as  well  the  difference  between  the  present  approach,  the 
incompressible  flow  approach,  and  the  fully  compressible  flow  approach  is  in  the 
assumption  regarding  the  equation  of  state.  The  compressibility  characteristics 
represented  by  the  equation  of  state  directly  affect  the  property  of  the  equation  of 
continuity  in  much  the  same  way  the  assumption  on  viscosity  affects  the  equation  of 
motion.  The  concept  of  compressibility  boundary  layer  greatly  speeds  up  the  computation 
of  small  Mach  number  flows  making  it  feasible  to  solve  complex  problem  with  currently 
available  computers.  The  large  eddy  simulation  concept  is  also  incorporated  here  in  order 
to  efficiently  simulate  the  dynamics  of  large  scale  eddies.  Specific  conclusions  concerning 
non-cavitating  flows  and  their  relationships  with  cavitating  flows  are  listed  below. 

1 .  When  the  cavitation  number  and  the  Reynolds  numbers  are  both  large,  the  non- 
cavitating  flow  pattern  depends  greatly  with  the  attack  angle.  For  small  angle  of  attack 
there  is  likely  to  be  no  boundary  layer  separation,  or  the  boundary  layer  separates  only  at 
the  trailing  edge.  The  separation  point  is  likely  to  move  upstream  as  the  angle  of  attack  is 
increased. 

2.  At  modest  angle  of  attack,  the  boundary  layer  may  breakup  and  produce  large- 
scale  vortex  shedding.  This  shedding  is  quite  periodic  if  the  foil  is  thin  and  the  separation 
point  is  fixed  at  the  leading  edge.  At  a  larger  attack  angle,  two  or  more  shed  eddies  may 
merge  near  the  trailing  edge  of  the  foil  producing  a  larger  and  stronger  eddy.  This  eddy 
washes  out  more  or  less  periodically  at  substantially  lower  fi’equencies  than  vortex 
shedding. 

3.  Lift  and  drag  oscillations  due  to  vortex  washout  and  the  wake  instability  are 
quite  significant.  In  contrast,  vortex  shedding  generates  relatively  weak  oscillations.  This 
is  due  to  the  fact  that  there  are  multiple  eddies  distributed  on  the  foil  surface,  and  their 
effect  on  forces  cancel  each  other  out. 

4.  The  vortex  shedding  and  washout  phenomena  of  hydrofoil  are  very  much 
similar  to  the  shedding  and  washout  phenomena  of  the  two-dimensional  straight  diffuser. 

5.  For  a  thin  foil  like  the  E.N.,  cavitation  is  likely  to  start  at  the  leading  edge.  The 
resulting  cavitation  is  quite  periodical  and  predictable.  Vorticity  generated  at  the 
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separation  point  continuously  feeds  the  sheet  cavity  to  make  it  grow  until  a  large  portion 
of  it  breaks  away  to  periodically  form  a  cloud  cavity.  As  the  cloud  cavity  moves  to  a 
higher  pressure  region  near  the  trailing  edge,  it  collapse  to  generate  very  high  pressure. 

6.  Detailed  computations  were  carried  out  for  the  case  of  a  6.2  degree  attack 
angle  and  the  cavitation  number  of  1.2,  and  the  results  were  compared  with  the 
experimental  results  of  Kubota  et  al.  The  computed  cavity  flow  patterns  qualitatively 
agree  well  with  the  experiment.  The  calculated  frequency  and  amplitude  of  pressure 
oscillation  on  the  foil  surface  due  to  the  dynamics  of  cloud  cavitation  agree  well  with  the 
experimental  data. 

7.  The  growth  and  breakdown  of  sheet  cavity  at  a  6.2  degree  attack  angle  are  very 
similar  to  the  leading  edge  flow  separation  and  vortex  shedding  of  non-cavitating  foil  at  a 
15  degree  attack  angle.  There  is  no  such  flow  separation  and  vortex  shedding  phenomena 
for  non-cavitating  flow  at  a  6.2  degree  attack  angle. 

8.  For  a  thick  foil  like  NACA0015  foil,  the  cavitation  pattern  is  not  as  regular  as 
that  of  a  thin  foil.  For  a  small  attack  angle  at  slightly  below  the  critical  cavitation  number, 
a  cavity  first  appears  near  mid  chord  which  travel  downstream  while  it  expands  and 
collapse.  This  cavitation  induces  boundary  layer  separation  and  vortex  shedding. 

9.  A  fixed  sheet  cavity  can  occur  near  the  nose  if  the  attack  angle  and  the 
cavitation  number  are  both  small.  No  cloud  cavity  is  produced  in  this  case. 

10.  For  a  modest  attack  angle  case,  two  flow  regimes  can  occur  intermittently. 
They  are:  a)  a  fixed  sheet  cavity  near  the  nose  and  b)  sheet  cavity  -  cloud  cavity  cycle. 

1 1 .  For  a  large  attack  angle,  it  is  most  likely  to  have  the  sheet  cavity  -  cloud  cavity 
flow  regime. 

12.  The  production  and  collapse  of  a  cloud  cavity  is  a  highly  3-dimensional 
phenomenon.  It  can  take  place  at  different  spanwise  locations. 
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Figures  4-1  through  6-23 
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Fig.  4-1  Medium  mesh  system  (381x64x3)  around  a  NACA0015  foil. 
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Attack=0.0  Deg,  Computed 
Attack=4.0  Deg,  Computed 
Attack=8.0  Deg,  Computed 


25 


Fig.  4-2  Comparison  of  computed  pressure  coefficients  with  experiment 
data  (Izumida,  et  al.,  1980)  for  NACA0015  foil  at  non-cavitating  condition. 


Fig.  4-3  Computed  pressure  coefficient  distributions  along  foil  surface 
based  on  three  different  mesh  systems:  coarse  mesh  (191x32x3),  medium  mesh 
(381x64x3),  and  fine  mesh  (761x128x3). 
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Fig.  4-4  Instantaneous  vorticity  component  ( ©r)  contours  on  the  mid-span 
cross-section  of  NACA662-415  foil  at  ct=6°  during  one  period  of  trailing  edge 
vortex  shedding. 
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Fig.  4-5  Velocity  (u)  spectra  at  5  positions  and  the  lift  spectra  for 
NACA662-415  foil  at  a=6“. 
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Fig.  4-7  Close-up  of  the  leading  edge  vortex  shedding  during  one  cycle 

featured  by  instantaneous  vorticity  component  contours  on  the  mid-soan 
cross-section  of  NACA662-4 1 5  foil  at  a= 1 0° . 
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Fig.  4-8  Velocity  (u)  spectra  at  5  positions  and  the  lift  spectra  for 
NACA662-415foUat  a=10°. 
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Power  Spectra 
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Fig.  4-10  Velocity  (u)  spectra  at  5  positions  and  the  lift  spectra  for 
NACA662-415  foil  at  a=15”. 


Fig.  4-11  Instantaneous  velocity  vectors  on  the  mid-span  cross-section  of 
NACA662-415  foil  at  a=l5\ 
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Power  Spectra 
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Velocity  (u)  Frequency  Analysis  by  FFT,  Attack=20  degrees 


Fig.  4-13  Velocity  (u)  spectra  at  5  positions  and  the  lift  spectra  for 
NACA662-415  foil  at  a=20°. 
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Fig.  5-2  Time-averaged  velocity  field  around  an  E.N.  foil  at  ^=6.2°  under 
non-cavitating  condition. 
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SHEET  AND  CLOUD  CAVITATION  SIMULATION 

Mean  Pressure  Coefficient  (Cp)  on  Foil  Surface  (EN  Foil),  Cl=0.754 


Fig.  5-5  Time-averaged  pressure  coefficient  distribution  along  E.N.  foil 
surface  at  a=6.2°  under  non-cavitating  condition. 


Unsteady  Cloud  Cavitation  (E.N.  Foil,  cc  =  6.2% 
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Fig.  5-6  Instantaneous  vorticity  component  (®z)  Fig.  5-7  Experimental  photographs  taken  by  Kubota 

contours  on  the  mid-span  cross-section  of  an  E.N.  foil  at  ot  al.  (1989)  for  an  E.N.  foil  at  oc=6.2  under  cavitating 

a-(i.T  under  cavitating  condition  of  <t=1  .2.  condition  of  <t=  1 .2. 


Velocity  Field  Comparison  (E.N.  Foil,  a  =  6.2% 


(b)  Experiment 

Fig.  5-8  Comparison  of  computed  velocity  field  with  experimental  result 
(Kubota  et  al.,  1989)  of  flow  over  an  E.N.  foil  at  a=6.2°  under  cavitating 
condition  of  tT=  1.2. 


Fig.  5-9  Instantaneous  velocity  field  around  an  E.N.  foil  at  a=6.2°  under 
cavitating  condition  of  0"=  1 .2. 
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Fig.  5-10  Close-up  of  instantaneous  velocity  fields  at  the  leading  edge  of  an 
E.N.  foil  at  (X=6.2°  under  cavitating  condition  of  cr=1.2,  featured  by  one  period 
leading  edge  vortex  shedding  due  to  the  cavitation. 


Fig.  5-11  Computed  time-space  (t,y)  distributions  of  vorticity  component 

(«z)  at  three  veitical  sections:  (a)  x/c=0.25;  (b)  x/c=0.65,  (c)  x/c=1.05,  for  an 
b.N.  toil  at  a=6.2  under  cavitating  condition  of  CT=1 .2. 


Fig.  5-12  Measured  time-space  (t,y)  distributions  of  vorticity  component 

( ®z)  at  x/c=1.05  (Kubota  et  al.,  1989)  for  an  E.N.  foil  at  a=6.2°  under  cavitating 
condition  of  cr=1.2. 
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(b)  Experiment 


Fig.  5-13  Comparison  of  computed  pressure  with  experimental  data 
(Kubota  et  al.,  1989)  at  three  locations  on  foil  surface  for  an  E.N.  foil  at  a=6.2° 
under  cavitating  condition  of  (T=1.2. 
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Pressure  Coefficient  (plus  Cavitation  No.) 


Fig.  5-14  Computed  1 
on  foil  surface  for  an  E.N.  foi 


iod  pressure  change  with  time  at  five  locations 
6.2°  under  cavitating  condition  of  cr=  1 .2. 
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v/r^-n  Companson  of  computed  nonnal  velocity  at  x/c=1.25  and 

y/^0. 133  with  expenmental  data  (Kubota  et  al.,  1989)  for  an  E.N  foil  at  a-6  T 
under  cavitatmg  condition  of  <t=  1.2. 
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SHEET  AND  CLOUD  CAVITATION  SIMULATION 

Mean  Pressure  Coefficient  (Cp)  on  E.N.  Foil  Surface,  Attack=6.2  degrees 


Fig.  5-16  Computed  time-averaged  pressure  distributions  on  E.N.  foil 
surface  at  a =6.2°  with  different  cavitation  numbers. 
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Pressure  Coefficient 


Fig.  6- 1  Relationship  between  pressure  and  density 
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(a)  High  speed  photograph  of  cavitating  flow 
taken  by  Kubota  et  al.  (1992). 


(b)  Computed  instantaneous  density 
distribution 


Fig.  6-2  Comparison  of  vortex-type  cavity  appearance  on  NACA0015  foil 
at  a=20°  under  cavitating  condition  of  <t=1.7,  (a)  experimental  photograph;  (b) 
computed  density  distribution. 
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Fig.  6-3 
ler  cavitating  ci 
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ressure  field  around  NACAi 
,0  (stable  sheet  cavity). 
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VORTICITY(ATTACK  ANGLE=8  DEG,  NACA0015) 
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Fig.  6-12  Instantaneous  pressure  field  around  NACA0015  foil  at  a=! 
under  cavitating  condition  of  (T=1.0  (sheet  and  cloud  cavities). 
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Fig.  6-16  Instantaneous  density  distribution  aroupd  NACA0015  foil  at 
=20°  under  cavitating  condition  of  (T=1.0  (sheet  and  cloud  cavities). 
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PRESSURE  FIELD  (ATTACK  ANGLE=8  DEG,  NACA0015.  NS.  NZ=29) 


Fig.  6-18  Instantaneous  pressure  field  around  NACA0015  foil  at  three 
span- wise  sections  at  a=8°  under  cavitating  condition  of  o-=1.0,  (a)  left  end 
section  (Section  1,  NZ=2),  (b)  mid-section  (Section  2,  NZ=15),  (c)  right  end 
section  (Section  3,  N2^28). 
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Section  1  (NZ=2) 


Section  2  (NZ=15) 


Section  3  (NZ=28) 


a=1.0  t/T=16.72 

DENSITY  FIELD  (ATTACK  ANGLE=8  DEG.  NACA0015.  NS.  NZ=29) 


Fig.  6-19  Instantaneous  density  distribution  around  NACA0015  foil  at 
three  span-wise  sections  at  a=8°  under  cavitating  condition  of  <t=1.0,  (a)  left  end 
section  (Section  1,  NZ=2),  (b)  mid-section  (Section  2,  NZ=15).  (c)  right  end 
section  (Section  3.  NZ=28). 
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span-wise  sections  at  a=8°  under 
section  (Section  1,  NZ=2),  (b)  mid-seci 
section  (Section  3,  NZ=28). 


DEG.  NACA0015.  NS.  NZ 


NACA0015  foil  at  three 
of  <T=1.0,  (a)  left  end 
NZ=15),  (c)  right  end 


Section  1  (NY=1,  Surface) 


Section  2  (NY=9) 


PRESSURE  FIELD  (ATTACK  ANGLE=8  DEG,  NACA0015,  NS.  NZ=29) 


Fig.  6-21  Instantaneous  pressure  field  around  NACA0015  foil  at  three 
vertical  levels  at  a=8°  under  cavitating  condition  of  (T=1.0,  (a)  foil  surface  level 
(Section  1,  NY=1),  (b)  y/c=0.05  (Section  2,  NY=9),  (c)  y/c=0.1  (Section  3, 
NY=15). 
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Section  1  (NY=1,  Surface) 


CT=1.0  1/T=16.72 

DENSITY  FIELD  (ATTACK  ANGLE=8  DEG,  NACA0015,  NS,  NZ=29) 


Fig.  6-22  Instantaneous  density  distribution  around  NACA0015  foil  at  three 
vertical  levels  at  a=8°  under  cavitating  condition  of  <t=1.0,  (a)  foil  surface  level 
(Section  1,  NY=1),  (b)  y/c=0.05  (Section  2,  NY=9),  (c)  y/c=0.1  (Section  3, 
NY=15). 
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Section  1  (NY=1,  Surface) 


(J=1.0  1/T= 16.72 

VELOCITY  FIELD  (ATTACK  ANGLE=8  DEG,  NACA0015.  NS,  NZ=29) 


Fig,  6-23  Instantaneous  velocity  field  around  NACA0015  foil  at  three 
vertical  levels  at  a=8°  under  cavitating  condition  of  (7= 1.0,  (a)  foil  surface  level 
(Section  1,  NY=1),  (b)  y/c=0.05  (Section  2,  NY=9),  (c)  y/c=0.1  (Section  3, 
NY=15). 
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